Dark soliton creation in Bose-Einstein condensates 
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It is demonstrated that stable, standing dark solitons can be created in current dilute-gas Bose- 
Einstein condensate experiments by the proper combination of phase and density engineering. Other 
combinations result in a widely controllable range of grey solitons. The phonon contribution is small 
and is calculated precisely. The ensuing dynamics should be observable in situ, i.e. without ballistic 
expansion of the condensate. 



The observation of solitons in Bose-Einstein conden- 
sates (BEC's) brings together two heretofore disparate 
fields of study: quantum degenerate gases [Q and ex- 
actly integrable nonlinear systems Q. Recently three- 
dimensional, solitonic structures have been observed 
in weakly-interacting atomic gas BEC's for the first 
time HQ. As fundamentally nonlinear, collective excita- 
tions of a macroscopic quantum wave they are intriguing 
in their own right. However, standing solitons in one di- 
mension, which are expected to have a wealth of special 
properties have not yet been created. 

The Gross-Pitaevskii equation which describes the 
mean field dynamics of a dilute BEC at low tem- 
peratures reduces to a one-dimensional nonlinear 
Schrodinger equation (NLS) when the transverse dimen- 
sions of the condensate are on the order of its healing 
length and its longitudinal dimension is much longer than 
its transverse ones ||. This is termed the quasi-one- 
dimensional (quasi- ID) regime of the Gross-Pitaevskii 
equation. If the transverse dimensions of the condensate 
are much less than the healing length then the Gross- 
Pitaevskii equation no longer applies and other physics 
ensue Q. The recent confinement of a BEC in a hollow 
blue-detuned laser beam demonstrates that the quasi- ID 
regime is experimentally realizable |^]. In three dimen- 
sions solitonic structures are unstable to transverse mod- 
ulations but in quasi- ID they are stable, and are in fact 
the stationary states of the NLS ]To[ . 

The extent to which quantum fluctuations affect the 
BEC is an outstanding theoretical question . As dark 
solitons are solutions to the mean field theory, substantial 
instability of dark solitons at low temperature in quasi- 
1D is an experimental measure of such higher order ef- 
fects. The present experimental technique for creating 
solitonic structures is based on phase engineering alone, 
makes additional transient density waves, and cannot be 
used to make a single standing dark soliton. Observa- 
tions have been made in three-dimensional harmonic ge- 
ometries or at temperatures for which about 10% of the 
atoms remained uncondensed. Time scales for observa- 



tion are subsequently on the order of 10 ms, whereas 
the lifetime of the BEC is on the order of seconds, pre- 
cluding the likelihood of observing the effect of quantum 
fluctuations. Finally, the length scale of these solitonic 
structures, i.e. the healing length, is so small that in 
situ observation of their dynamics is not possible and 
the condensate must first be expanded by turning off the 
confining potential. In the following, a method which 
remedies all of these difficulties is presented. 
The NLS may be written ]lo| in the form 

[~ed xx + L\xP(x,t)\ 2 + V(x)}xP(x,t) =idt1>(x,t), (1) 

where ip is the macroscopic quantum wavefunction, 
L is a longitudinal confinement length, V(x) = 
(2m^ 2 /h 2 )V (x) is a rescaling of a confining potential 
V (x), £ = (87rna) -1 / 2 is the healing length, a is the 
s-wave scattering length for binary atomic interactions, 
and n is the mean number density, ip is normalized to 
1 and has units of L -1 / 2 , x has units of length, and 
t = (2m£ 2 /h)t Q is dimensionless, where t Q has units of 
time. 

In particular, we consider the case of 87 Rb confined in 
a box-like potential, with L = 100 jam and a transverse 
length of L t = 10 /im. For N = 1.1 x 10 4 atoms and 
n = N/(LL 2 ) the effective healing length [[l2| is 2.5 /im. 
As the length scale of solitons is ~ 2£ and the wave- 
length of the imaging radiation is ~ 0.5 /im, this exper- 
imentally realizable configuration ensures observability 
without needing to ballistically expand the condensate. 
It also ensures the possibility of phase and density engi- 
neering of structures on the scale of the healing length 
by means of the dipole potential of laser fields. A key 
point in this choice of parameters is the use of a box-like, 
rather than harmonic, confining potential For har- 
monic confinement the Thomas-Fermi radius jlj of the 
condensate scales as A^ 1 / 5 , n oc N/N 3 / 5 , and therefore 
£ oc A/ --1 / 5 , whereas for box-like confinement £ oc A" -1 / 2 . 

Neglecting the effect of V(x) and writing the wavefunc- 
tion as ip(x,t) = y/ p(x, t) exp[i0(x, t) — ifit], the single 
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FIG. 1. (a) A combination of a box-like potential and a 
tightly focused, blue-detuned laser beam is used to engi- 
neer the density, (b) The resulting wavefunction is phase 
engineered with a second, far-detuned laser beam, resulting 
in an initial state very close to a standing dark soliton. 



soliton solution to Eq. ([!]) takes the form 

p(x,t) = - {l-2j 2 £ 2 sech 2 [j(x- ct)}} (2a) 



4>{x,t) 

where 



1 | ^ tanhfr (x - c t)} j + | (2b) 



tan 

2 7 2 £ 2 + c 2 /(2£ 2 ) = 1 . 



(2c) 



For L ^> £ Eqs. (g) are an excellent approximation away 
from the walls. The single particle density p(x,i) is 
constant except over a density notch of width 1/7, and 
the phase changes sharply and monotonically across the 
notch. The chemical potential fi = 1 in these units 
and c is the speed of the soliton. The constraint (||c) 
links the soliton velocity c to its depth 2£ 2 7 2 and thus 
leaves only a single free parameter for the soliton solu- 
tion. In the case c = 0, I/7 = ^/2( i and the soliton 
forms a node, is a stationary solution to Eq. ([!]), and 
is called a dark soliton. In the case < c < c s , where 
c s = y/2£ is the Bogoliubov sound speed, the soliton is 
a moving density notch and is called a grey soliton. In 
the case c — > c~ , 1/j diverges and the soliton depth ap- 
proaches zero. Upon reinsertion of constants it is found 
that c s = V47rna %/m. The total phase difference across 
the notch is A0 = 2tan- 1 [ v ' / l - (c/c s ) 2 / (c/c s )}. 

As is apparent from Eq. (||), a soliton requires both 
a nontrivial phase and a nontrivial density profile. The 
technique of phase engineering used to make solitonic 
structures in previous experiments |||| imprinted a con- 
stant phase on one half of a nearly uniform condensate. 
This resulted in a combination of transient density waves 
and one or more grey solitons. We propose instead the 
following scheme. 

Atoms are condensed into a box-like trap with a 
sharply focused, blue-detuned laser beam in the cen- 
ter of the long trap axis. Due to a large laser detuning 
from atomic resonances, spontaneous processes can be 
neglected on the relevant time scales. Thus the effect of 
the laser beam on the condensate atoms can be described 
by the optical dipole potential, which is proportional to 
the local light intensity. In Fig. |l|(a) is shown the poten- 
tial used as a model, a gaussian with a width on the order 
of the healing length, i.e. 2Ax v ~ 2.5/im, steep walls at 
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FIG. 2. A standing dark soliton is created with the 
proper combination of phase and density engineering. In 
addition to the central dark soliton, characterized by a 
7r phase difference across a density notch which forms a 
node, a spray of phonons is emitted at the speed of sound 
(c s = a/2) to either side. A second, shallow soliton, emitted 
to the left at nearly c s , offsets the effect of the finite slope 
of the initial phase due to the diffraction-limited fall-off of 
the phase-engineering laser beam. 



±L/2 , and a height on the order of the chemical poten- 
tial fi = 1 : 
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(3) 



In Fig. [l] and in all simulations the parameters £ = 1 
and L = 40 were used. Note that it should also be pos- 
sible to create a density notch by adiabatically ramping 
the intensity of a laser into an initially uniform BEC. 
The resulting deformation dips to ~1% of the maximum 
density, or ~10% of the maximum amplitude, as may 
be seen in Fig. |l|(b). Because the scale of variation of 
V(x) is < £ the kinetic energy term in Eq. (0) cannot 
be neglected and the response of the condensate is not in 
the Thomas- Fermi regime. It is therefore obtained nu- 
merically by imaginary time relaxation of Eq. (||) with 
Eq. ®. 

At t = the focused laser beam is switched off and a 
second, far-detuned laser pulse of uniform intensity dis- 
tribution is shined on one half of the condensate. The 
pulse duration is chosen to be shorter than the correla- 
tion time of the condensate, t c = £/c s . This ensures that 
the light field changes only the phase distribution and not 
the density distribution of the BEC. In the simulation it 
therefore suffices to switch on a controllable phase (j){x) 
instantaneously M. 



4>{x) = tanh(2x/Ax0) 



(4) 



was used as a model, which represents the diffraction- 
limited fall-off of the laser beam over the notch, with a 
width of and a total phase difference of A0, as may 
be seen in Fig. |l|(b), where the parameters A</> = tt and 
Ax^ = £ ~ 2.5 /im have been used. 

The resulting notch is both density and phase engi- 
neered. In Fig. is shown its evolution. For 87 Rb there 



2 




-20 -10 10 



-20 -10 10 x 



FIG. 3. A variety of physically intriguing possibilities 
arise from variations in the technnique used to create a 
standing dark soliton. (a)-(b) The use of density engineer- 
ing alone creates a pair of equal and opposite grey solitons 
with a widely controllable range of velocites. (c)-(d) Phase 
engineering with a phase different from ty creates an asym- 
metric grey soliton pair. Note the spatial shift inherent in 
soliton-soliton interactions. 

are 2.27 ms per simulation time unit ]To| |. When Acf) = tt 
was used it was found that the soliton drifted slowly to 
the left, i.e. was not completely dark. A true stand- 
ing dark soliton requires a step function in the phase. 
Therefore to counter the effect of the diffraction-limited 
fall-off of the second laser beam, A<j) = 1.05 x tt was used. 
This creates a second, shallow grey soliton which carries 
away the small momentum of the first without apprecia- 
bly deforming the density. Thus a standing dark soliton 
is successfully created. 

As is apparent in Fig. ||, some portion of the initial 
notch is radiated as phonons. The apportionment of the 
total energy into phonons and solitons may be calculated 
precisely. Because the NLS is nonlinear it is not possi- 
ble to build any particular solution by a summation over 
phonon modes. As phonons and solitons are the funda- 
mental solution types of the ID NLS it is therefore use- 
ful to single out their contribution to the dynamics. The 
NLS has a denumerably infinite number of integrals of 
motion || . In the context of the BEC the first three cor- 
respond to normalization, energy, and momentum. For 
the simulations studied herein it suffices to consider en- 
ergy. 



The second integral of motion of Eq. (|l|) is 



,L/2 

/ dx 
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e\d x ^\ 2 + -m 4 



V[x) 



Assuming L ^> £, the energy of a single soliton is 



E s ~N 



dx 



(5) 



(6) 



where Eq. (||) has been renormalized to subtract out the 
constant background and the potential V(x) has been as- 
sumed to be constant in the region of the soliton. The 
factor of N follows from the use of a single, rather than 
TV, particle density for \ip(x, t)\ 2 . Equations (|^) may then 
be used to obtain the energy of a single soliton: 



N 4 
L3 1 



|sin(A0/2)| 2 



(7) 



E s is dimensionless and depends on the line density N/L. 
The rescaling [h 2 / (2m£ > 2 )]E s oc (N/L) 3 / 2 reinstates the 
units of energy. For our choice of parameters the differ- 
ence between the approximation and an exact expression 
for E s is negligible and not observable in simulations. 

Conservation of the total energy of the system may be 
expressed as Ed + E Q = E s + E p + E Ql where E&, E Ql 
and E p are the energies of the initial deformation, the 
undeformed background ground state, and the phonons, 
respectively. By use of Eq. (||) Ed and E Q may be de- 
termined in the simulation. E s may be calculated from 
Eq. (0). With this method it was determined that the 
total phonon contribution in Fig. || is 39% of the defor- 
mation energy Ed and the second soliton carries away 
only 1.5% of Ed- However, the phonons spread out over 
the whole box while the 60% of the energy carried by the 
dark soliton is localized in a region of width ~ 2£ = 2. 
Therefore this method of dark soliton creation is highly 
efficient and phonons are essentially small fluctuations in 
the density. 

It is important to explore the robustness of the pro- 
posed technique. To what extent does it depend on the 
parameters of Eqs. (||)-([f)? To this end many simulations 
were performed which varied Ax v , Ax^, and A<p. It was 
found that 2Ax v < £ was optimal. Note that the full 
width of the gaussian is 2Ax v . 2Ax v = 2£ produced an 
additional pair of very shallow grey solitons which would 
not be observable in experiments. Ax^ < £ was sufficient 
to produce a standing dark soliton while L ^> Ax^ > £ 
set the soliton in motion. Finally Acf) ~ tt was required, 
as expected. It was also found that the amount of the 
deformation energy transformed into phonons was con- 
sistently about 40% to 50%. Outside of these parame- 
ter ranges a number of additional intriguing possibilities 
arose. In the following a few such possibilities are illus- 
trated. 

If density engineering alone is used, a pair of equal and 
opposite grey solitons emerge. This effect is illustrated in 
Fig. ||(a)-(b). The velocities of the grey soliton pair may, 
by use of Eq. (^), be used to determine the total energy 
of the phonons. It was found in simulations that within 
the range of 2Ax = 0.5 to 2.5 from 40% to 70% of the 
deformation energy was converted into phonons and the 
remaining energy was converted into a symmetric grey 
soliton pair. In general a lower energy notch resulted 
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FIG. 4. (a)-(b) Density engineering a notch wider than 
the healing length creates a symmetric fan of grey solitons. 
(c)-(d) If phase engineering is added a dark soliton is cre- 
ated in the center of the fan. 



noted that solitons are robustly stable structures (T(| 
and, though not shown here, the introduction of 10% 
stochastic noise into the simulations had no observable 
effect on the results. 

In conclusion, a method of creating a single, standing 
dark soliton in a gaseous Bose-Einstein condensate has 
been presented. This experimentally realizable method 
combines phase and density engineering and requires one 
or two lasers and a quasi-one-dimensional trap. Varia- 
tions result in the creation of grey solitons with a widely 
controllable range of velocities. 

We thank Klaus Sengstock and William Reinhardt for 
helpful discussions. This work was supported by National 
Science Foundation Grant CHE97-32919 and by SFB407 
of the Deutsche Forschungsgemeinschaft. J. Brand ac- 
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in slower solitons, with a lower limit on the soliton ve- 
locity of ~ 0.45c s . In previous solitonic-structure ex- 
periments, which used phase engineering alone ||, the 
velocities were 0.54c s to 0.76c s and 0.68c s to c s , respec- 
tively HQ . By the method of density engineering a wide 
range of velocities from 0.45c s to c s may be obtained by 
adjusting the width or depth and thereby the energy of 
the initial deformation. 

Another variation on the technique, which uses both 
density and phase engineering and a phase difference of 
Ac/) 7^ 7r, creates an asymmetric grey soliton pair. 60% 
to 70% of the deformation energy is radiated as phonons. 
For the phase profile of Eq. (|^), 7r < Acj) < 2n results in 
a faster soliton to the right. < Acj) < tt results in a 
faster soliton to the left, as is to be expected from the 
constant phase offset inherent in the wavefunction. The 
total phase difference across the two solitons adds up to 
A(j). Such asymmetric interactions are useful in observing 
the delay inherent in soliton interactions ]l5[ . 

If again density engineering alone is used, but with too 
large a notch, i.e. Ax v > £, multiple pairs of equal and 
opposite grey solitons emerge, as has also been predicted 
elsewhere |14|| . If phase engineering is added to this sce- 
nario a central, dark soliton is created in the middle of 
the trap with equal and opposite pairs of grey solitons 
fanning outwards to either side. These cases are shown 
in Fig. f| 

In general the consistent 40%- 70% of the deformation 
energy radiated into phonons is representative of the dif- 
ference in form of the density deformation which results 
from the gaussian in Eq. (||) and the sech in Eq. (|2|a). 
A different density notch could be expected to result in 
a different apportionment of energy. It should also be 
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